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ABSTRACT 

We consider a forming galaxy undergoing multitudinous supernova (SN) explosions, as a possible 
model of Lya blobs (LABs). For this purpose, an ultra-high resolution hydrodynamic simulation is 
performed using 1024^ grid points, where SN remnants are resolved with sufficient accuracy. It is found 
that multiple SN explosions produce kpc-size expanding hot bubbles, which drive cool, dense shells by 
strong shock. The colliding high-density cooling shells radiate intensive Lya emission, resulting in a high 
Lya luminosity of ~ lO'''^ erg s^^, comparable to the observed level in LABs. Also, recently discovered 
bubbly features in some LABs are quite similar to the structure predicted in the present simulation. 
Furthermore, the result demonstrates that LABs are representative of evolving primordial galaxies; they 
could hold direct information on the early chemical enrichment of galaxies, contrary to present-day 
galaxies which have undergone intense recycling of interstellar matter, thus erasing most of the early 
chemical history. It turns out that the metal mixing proceeds in a very inhomogeneous fashion, so that 
there appears a large spread of metallicity, that is,[Fe/H] « to — 5 or [O/H] w 1 to — 4. Hence, the early 
galactic chemical evolution may have proceeded in a different manner from that hitherto considered in 
one-zone models. 

Subject headings: cosmology: theory — galaxies: formation — galaxies: high-redshift — galaxies: 
starburst — hydrodynamics — supernovae: general 



1. INTRODUCTION 

Lya blobs (LABs) are a new category of high redshift 
galaxies, which have been recently discovered. Steidel et 
al. (2000) found two bright Lya blobs (Blob 1 and Blob 
2) in the SSA22 field, which have sizes of ~ 100 kpc and 
Lya luminosities of ~ lO^'^ erg s"""^ at z — 3.09. These 
LABs are associated with previously known Lyman-break 
galaxies (LBGs) at redshift (z) ^ 3.09. Compared to typ- 
ical Lya emitters found at high redshifts (Dey et al. 1998; 
Weyman et al. 1998; Rhoads et al. 2000; Hu et al. 2002; 
Ajiki et al. 2002; Dawson et al. 2002; Fujita et al. 2003; 
Kodaira et al. 2003; Matsuda et al. 2004), LABs are 
more luminous and much more extended. LABs resemble 
the giant Lya nebulae associated with high-redshift radio 
galaxies, although the associated radio continuum flux is 
less than 1%. Recently, Matsuda et al. (2004) have found 
35 candidates for LABs in the SSA22 field. One third of 
them are apparently not associated with UV continuum 
sources that are bright enough to produce Lya emission. 
Therefore, the origin of LABs is mysterious. 

So far, three possibilities have been proposed as models 
for LABs. The first is heavily-obscured UV sources like an 
active galactic nucleus (AGN) or starburst (Steidel et al. 
2000; Chapman et al. 2004). Chapman et al. (2004) ar- 
gue that the multiwavelength observations of Blob 1 is at 
least consistent with a dust-enshrouded AGN surrounded 
by a starburst. The second is a cooling radiation model 
(Haiman, Spaans, & Quataert 2000; Fardal et al. 2001). 
In a bottom-up scenario of galaxy formation, a large num- 
ber of subgalactic halos are expected to collapse at high 
redshifts and they can emit significant Lya fluxes through 
collisional excitation of hydrogen. Thus, the cooling ra- 



diation from proto-galaxies have been advocated to be a 
possible mechanism for observed Lya emission of LABs. 
The last is a superwind model (Taniguchi & Shioya 2000; 
Ohyama et al. 2003). Taniguchi & Shioya (2000) have 
argued that a large-scale bipolar outflow driven by the 
starburst in a forming galaxy can be responsible for the 
intensive Lya emission of LABs. 

The deep optical spectroscopy of LABs suggests that 
the kinematical properties of LABs favor supernova (SN) 
driven winds (Ohyama et al. 2003). Interestingly, Mat- 
suda et al. (2004) have recently revealed bubbly features in 
Blobs 1 and 2, where the bubble size is typically ~ 15kpc. 
The bubbly structures strongly suggest that SN events 
could be closely related to LABs. Therefore, complexes 
of various superbubbles, which can be driven by multiple 
SNe, are an attractive explanation for LABs. However, it 
has not been shown whether multiple SN explosions can 
actually produce the observed LAB emission. In this Let- 
ter, we attempt to build up a multiple SN explosion model 
for LABs. For the purpose, we perform an ultra-high res- 
olution (1024^^ grids) simulation on the SN explosions in 
inhomogeneous and multiphase interstellar matter (ISM), 
where individual SN remnants are resolved. The results of 
this chemodynamical simulation, incorporating spectro- 
photometric modeling (see Mori et al. 1997 and Mori, 
Yoshii & Nomoto 1999), can be directly compared to the 
observations on Lya blobs. In §2, the numerical model is 
described. In §3, we present the numerical results and in 
§4 the early chemical evolution of the system is discussed. 



2. MODEL 
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We define a fiducial model for a Lya blob at redshift 
2 = 3, in which the total mass of gaseous matter is ini- 
tially assumed to be Mg = 9.5 x 10^^ M©, and the total 
mass including dark matter halo is Mt = 3.7 x 10^^ Mq. 
The baryonic matter to dark matter halo has been in- 
creased by a factor « 2 with respect to the cosmological 
value*, to account for possible bias effects in rapidly cool- 
ing regions. According to the general picture of bottom- 
up scenarios for galaxy formation, we model a protogalaxy 
as an assemblage of mimcrous subgalactic condensations 
building up the total mass of a galaxy. As for the surface 
density profiles of the Lya blobs, the galaxy as a whole 
is modelled as a pre-virialized system with a spherically- 
averaged density profile pd{r) oc /Onfw('')^^^ in a radius 
of 106 kpc, where pnfw(^) is a NFW density profile, 
PNFw(r-) oc r"i(r + a)"2 (Navarro, Frenk & White 1997). 
Subgalactic condensations having a zero velocity are then 
distributed to such a profile by Monte-Carlo realization. 
These are taken as already virialized systems with a NFW 
density profile encompassing a mass of 1.4 x 10^ M© in a 
virial radius Vg = 8.60 kpc. 

The star formation rate (SFR) is assumed to be propor- 
tional to the total gas mass. The SFR integrated over 
the whole system is SFR(t) = 204exp(-t/t^) Moyr"^ 
with the time scale = 4.6Gyr. Because of the short 
timescale of our simulation, the SFR is essentially con- 
stant ~ 200 MQyr~^ which is inferred from observations 
of bright LBGs at z - 3 (Shaplcy et al. 2001). Thus, the 
ambiguity of has little serious problem. The mass of 
individual stars and the formation epoch are determined 
by randomly sampling both SFR(i) and the initial mass 
function (IMF), which is assumed to be of Salpeter type 
with upper and lower mass cutoffs of 120 and 0.1 Mq, 
respectively. The newly born stars are distributed in pro- 
portion to the gas density using a Monte-Carlo procedure. 
A star more massive than 8 Mq is assumed, after the main 
sequence lifetime, to undergo a Type II SN explosion, re- 
leasing a total energy of 10''^ erg. A mass of 2.4M0 of 
oxygen and of 9.05 x 10~^Mq of iron are ejected from 
a Type II SN explosion (Tsujimoto et al. 1995). In the 
present simulation, a total of 7.53 x 10^ SN explosions oc- 
cur in the whole galaxy; we focus on the first 50 Myr of 
the chemodynamical evolution of the system. 

The gas dynamics is solved in three-dimensional space 
with 1024^ Cartesian grids. Here, we employ the AUSM- 
DV scheme (Wada & Liou 1994), which is a TVD scheme 
that can resolve the shock accurately. The spatial res- 
olution is 195 pc and the simulation box has a (physi- 
cal) size of 212 kpc. The ejected energy and well-mixed 
metals by SN are supplied to 8 local cells surrounding 
SN region. Then, the hydrodynamic evolution of met- 
als is followed by the same algorithm as the gas den- 
sity. Since the time variation of dark matter halo po- 
tential is insignificant in the early stage of galaxy evo- 
lution, we assume the gravitational potential to be con- 
stant in time. Also, the gas is assumed to be optically 
thin and in coUisional-ionization equilibrium. Radiative 
cooling is included self-consistently with matallicity, using 
the metallicity-dependent cooling curves by Sutherland & 
Dopita (1993). The set of basic equations are numerically 

^ Throughout this paper we assume a ACDM cosmology with Qm = 0. 
and n^fi^ = 0.02. 



solved by a parallel version of the Astrophysical Fluid Dy- 
namics (AFD2) scheme (see Mori, Ferrara & Madau 2002). 

3. RESULTS 

Figure 1 shows the numerical results at 50 Myr. Figures 
la and Ic show the density and temperature distributions 
in a slice along the X — Y plane, respectively. After the 
first massive stars explode, the gas temperature is raised 
up to about 10^ K locally, and expanding hot bubbles of 
~kpc are produced. Then, the strong shock drives cooled, 
dense shells enclosing hot bubbles. Subsequent SN explo- 
sions further accelerate the expansion of hot bubbles and 
the ambient gas is continuously swept up by the dense 
shells; the gas density in shells increases further owing to 
the enhanced cooling. Since the density of ISM is low 
in outer regions, the bubble expansion is faster there and 
SN-driven shock waves quickly collide with each other to 
generate larger bubbles of ~10kpc, which are surrounded 
by high density ( ^ 0.1 cm~^), cool (^ 10"* K) shells. 

We compute the emission properties of the galaxy as- 
suming an optically thin gas in coUisional ionization equi- 
librium and using the MAPPINGS III code by Suther- 
land & Dopita (1993). Lya extinction due to dust seems 
negligible. For the average gas column density, Nh ~ 
5.6 X 10^^ cm~^, and metallicity in the simulation (shown 
below), we have E{B - V) = (7Vh/9.2 x lO^^cm'^) x 
;LO[Fo/h] ^ 3 Qg X ;^Q-3 rpj^g^^ ^gjj^g equation (20) of Pei 
(1992) and Ab = 1.2 x lO'^, we find TLya = 3.1 x lO^^. 
Figures le and If show the projected distribution of the re- 
sultant Lya emission and X-ray emission in the 0.2-2 keV 
band (0.8-8 keV at 2 = 3), where one arcsec corresponds to 
8kpc. It is worth noting that small bubbles with a few kpc 
are smeared out in the surface brightness of Lya emission 
and bubbly features of ~ lOkpc are prominent in outer re- 
gions. Such features are quite similar to the bubbly struc- 
tures observed in Blobs 1 and 2 (Matsuda et al 2004). The 
total flux of Lya emission is 1.2 x 10^^® erg cm^^ and 
the corresponding luminosity is 9.7 x 10^^ erg s~^, which 
comes mainly from high-density cooling shells. This Lya 
luminosity is also comparable to ~ 10*^ erg s~^ of Blobs 1 
and 2. These results suggest that LABs could be an early 
phase of chemical evolution in a protogalaxy. The total 
flux of X-ray emission is 1.5 x 10~^^ erg s~^ cm~^ and 
the corresponding luminosity is 1.1 x 10*" erg s~^. This is 
within the constraint of the X-ray luminosity from Chan- 
dra X-ray observation for Blob 1, Lx < 2.1 x 10*^ erg s~^ 
(Chapman et al. 2004). 

The line-of-sight velocity distribution is also shown in 
Figure Id. Since gas in outer regions has a velocity higher 
than 500 km s"-*^, which exceeds the escape velocity of this 
galaxy (472 km s^^), it may escape from the galaxy poten- 
tial well, while in inner regions gas velocity is lower than 
the escape velocity, and therefore gas is confined within 
the gravitational potential. If the total mass of galaxy is 
lower, winds will be more effective (Mori, Yoshii & Nomoto 
1999; Ferrara, Pettini & Shchekinov 2000), and provide an 
efficient mechanism to distribute heavy elements over cos- 
mological volumes (Mori, Ferrara & Madau 2002). In the 
future, the more detailed bubbly structure seen in Figures 
le may be directly detected by high resolution observa- 

.3, Ha = 0.7, a Hubble constant of /i = i?o/(100 km s'^ Mpc-^) = 0.7 
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tions (~ 0.1 arcsec), of redshifted emission (1.5 x 10^^ 
erg s~^), with the James Webb Space Telescope. 

4. DISCUSSION 

As shown in the previous section, LABs could corre- 
spond to a quite early phase of galaxy chemical evolution. 
In this case, the observed bubbly features imply that the 
self enrichment is on going in a very inliomogcneous fash- 
ion in LABs. Actually, in the present simulation, the dis- 
tribution of metallicity is very inhomogeneous as shown in 
Figure Ic, where the predicted oxygen abundance [0/H] 
is shown. The gas in the vicinity of SNe is polluted with 
newly synthesized heavy elements ejected from SNe, hut a 
large amount of the gas still retains of low metallicity. The 
interactions of expanding hot bubbles give rise to a com- 
plex structure in the inner regions, where a metal-rich gas 
[0/H] ^ — 1.74 coexists with an almost primordial gas. 
They are separated from each other by cool shells. Inter- 
estingly, the outer regions tend to have higher metallicity 
compared with inner regions. This is because the metal en- 
riched super bubbles arc less diluted with metal-deficient 
gas in the outer low-density regions. 

Figure 2 shows the relation between density and metal- 
licity, where the abundance is represented by [Fe/H] as 
well as [0/H] and the contours denote the mass level. 
This figure clearly shows that the mixing of heavy ele- 
ments proceeds incompletely in an very inhomogeneous 
manner. The gas with higher metallicity of [0/H] > 
([Fe/H] > —0.5) is distributed only in low density re- 
gions of rig ^ lcm~^, which corresponds to hot super- 
bubbles. On the other hand, a part of the SN ejecta mix 
with the dense regions with primordial abundance (either 
shells or overdensities around subgalactic condensations) 
and tend to be diluted to lower metallicity of [0/H] ^ — 1.7 
([Fe/H] ~ —2.3). Note that considerable dispersion of 
metallicity (—5 ^ [0/H] ^ 0) is found in high density re- 
gions of rig ^ Icm"'^. This spread of metallicity hardly 
depends upon the numerical resolution as long as hot bub- 
bles are fully resolved in the simulation. Hence, the next 
generation of stars which are expected to form out of the 
high-density regions should exhibit a similar dispersion of 
metallicity. 

So far, the theoretical models of galactic chemical evo- 
lution have often assumed the homogeneous ISM (one- 
zone model), with the instantaneous and perfect mixing 
of heavy elements synthesized in SNe. On the contrary, 
the present simulations show that the metal mixing pro- 
ceeds in a very inhomogeneous fashion. This means that 
one should carefully treat the metal mixing for the study 
of chemical evolution of galaxies. 

The volume-averaged metallicity of the gas in the 
present simulation is found to be [Fc/H]v = —2.30, while 
the density-weighted metallicity is [Fe/H]p = —2.51. On 



the other hand, if one derives the mean metallicity based 
on X-ray flux, [Fe/H]x = —0.71 is found, which is about 
forty times higher than the volume-averaged one. This 
is because X-ray is mainly emitted from the metal-rich 
gas in the hot bubbles. The density-weighted metallic- 
ity is well above a so-called critical metallicity, Zcrit ~ 
IQ-^^^Zq (Omukai 2000; Bromm et al. 2001; Schneider et 
al. 2002) marking the transition from a top-heavy IMF to 
a normal (i.e. Salpeter) IMF. Thus, the next stellar gener- 
ation will form with an IMF similar to the present day one, 
with the low-mass stars possibly forming the metal-poor 
halo population (McWilliam et al. 1995; Ryan, Norris & 
Beers 1996). In addition, it is worth noting the existence of 
very low metallicity ([Fc/H] « -5), high density (rig ~ 1 
cm~^) regions. Since the abmidancc of these regions can 
be below the critical metallicity, the IMF could be similar 
to that of Pop III stars (e.g. Nakamura & Umemura 2001). 
Recently, an extremely metal poor star ([Fc/H] « —5.3) 
is discovered by Christlieb et al. (2002) may be a relic 
of such stars. It should be kept in mind, however, that 
the model galaxy is constructed from primordial gas in 
this paper. At redshift z = 3, Lya absorption lines in 
quasar spectra show that at least some fraction of the 
IGM is polluted to 10~^-^ Zq. If a galaxy forms from a 
well-mixed pre-enriched medium, the metallicity would 
be accordingly shifted to higher values, hence precluding 
the formation of extremely metal-poor stars. 

5. CONCLUSIONS 

We have suggested that LABs can be identified with pri- 
mordial galaxies catched in a supernova-dominated phase. 
This conclusion is based on ultra-high resolution hydrody- 
namic calculations which resolve individual SN remnants. 
As a result, we found that the resultant Lya luminosity 
can account for the observed luminosity of LABs; in ad- 
dition, the bubbly structures produced by multiple SN 
explosions are quite similar to the observed features in 
Lya surface brightness of LABs. Hence the emerging the- 
oretical picture is fully consistent with experimental data. 
If this picture is correct, bubbly features may show that 
the metal mixing is highly inhomogeneous in LABs. The 
present simulation predicts a quite wide range of metallic- 
ity for the next generation of stars. This may affect the 
galaxy evolution significantly. 
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REFERENCES 



Ajiki, M. ct al. 2002, ApJ, 576, L25 

Bromm, V., Fcrrara, A., Coppi, P. S., & Larson, R. B. 2001, 

MNRAS, 328, 969 
Chapman, S. C, Scott, D., Windhorst, R. A., Prayer, D. T., Borys, 

C, Lewis, G. F., & Ivison, R. J. 2004, ApJ, 606, 85 
ChristUeb, N. et aL 2002, Nature, 419, 904 



Dawson, S., Spinrad, H., Stern, D., Dey, A., van Breugel, W., de 

Vrics, W., & Rculand, M. 2002, ApJ, 570, 92 
Dcy, A., Spinrad, H., Stern, D., Graham, J. R., & Chaffee, F. H. 

1998, ApJ, 498, L93 
Fardal, M. A., Katz, N., Gardner, J. P., Hernquist, L., Weinberg, 

D. H., & Dave, R. 2001, ApJ, 562, 605 
Ferrara, A., Pettini, M., & Shchekinov, Y. 2000, MNRAS, 319, 539 



4 



Mori, Umemura & Ferrara 



Fujita, S. S. et al. 2003, AJ, 125, 13 

Haiman, Z., Spaans, M., & Quataert, E. 2000, ApJ, 537, L5 

Hu, E. M., Cowie, L. L., McMahon, R. G., Capak, P., Iwamuro, F., 

Kneib, J.-R, Maihara, T., & Motohara, K. 2002, ApJ, 568, L75 
Kodaira, K. et al. 2003, RASJ, 55, L17 
Matsuda, Y., et al. 2004, AJ, 128, 569 

McWilliam, A., Rreston, G. W., Sneden, C, & Searle, L. 1995, AJ, 

109, 2757 

Mori, M., Ferrara, A., & Madau, P. 2002, ApJ, 571, 40 
Mori, M., Yoshii, Y., & Nomoto, K. 1999, ApJ, 511, 585 
Mori, M., Yoshii, Y., Tsujimoto, T., & Nomoto, K. 1997, ApJ, 478, 
L21 

Nakamura, F. & Umemura, M. 2001, ApJ, 548, 19 

Navarro, J. F., Frank, C. S., & White, S. D. M. 1997, ApJ, 490, 493 

Ohyama et al. 2003, ApJ, 591, L90. 

Omukai, K. 2000, ApJ, 534, 809 

Rei, Y. C. 1992, ApJ, 395, 130 



Rhoads, J. E., Malhotra, S.. Dcy, A., Stern. D.. Spinrad, H., & 

Jannuzi, B. T. 2000, ApJ, 545, 85 
Ryan, S. G., Norris, J. E., & Beers, T. C. 1996, ApJ, 471, 254 
Schneider, R., Ferrara, A., Natarajan, P., & Omukai, K. 2002, ApJ, 

571, 30 

Shapley, A. E., Steidel, C. C, Adelberger, K. L., Dickinson, M., 

Giavaliseo, M., & Pettini, M. 2001, ApJ, 562, 95 
Steidel, C. C, Adelberger, K. L., Shapley, A. E., Pettini, M., 

Dickinson, M., & Giavaliseo, M. 2000, ApJ, 532, 170 
Sutherland, R. S. & Dopita, M. A. 1993, ApJS, 88, 253 
Taniguchi, Y. & Shioya, Y. 2000, ApJ, 532, L13 
Tsujimoto, T., Nomoto, K., Yoshii, Y., Hoshimoto, M., Yanagida, 

S., & Thielemann F.-K. 1995, MNRAS, 277, 945 
Wada, Y. & Liou, M. S. 1994, AIAA Paper, 94-0083 
Weymann, R. J., Stern, D., Bunker, A., Spinrad, H., Chaifee, F. H., 

Thompson, R. I., & Storrie-Lombardi, L. J. 1998, ApJ, 505, L95 



The Nature of Lya Blobs 



5 



-i --i -S i * i » 1 B -SSi -» -Sii -tl 




Fig. 1. — Numerical results at 50 Myrs on the multiple SN explosions in a forming galaxy. Here are shown the sectional distributions of (a) 
the density, (b) the oxygen abundance [O/H], (c) the temperature, and (d) the line-of-sight velocity in a slice along the J^-F plane. Also, 
(e) the projected distribution of the Lya emission and (f) that of the soft X— ray emission in the 0.2-2 keV band (0.8—8 keV at the redshift 
z=3) are shown. Here the angular resolution is 2.7 X 10~^ arcsec (corresponding to the numerical resolution, 195 pc). 
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Fig. 2. — Relation between gas density and metallicity. The contours show mass levels in logarithmic scales in units of M©. Oxygen 
abundance [O/H] corresponds to iron abundance of [Fe/H] =[0/H] -0.559 for the yield adopted. 



